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The Josephson current flowing through weak links containing ferromagnetic elements is studied 
theoretically under the condition that the quasiparticle distribution over energy states in the in- 
terlayer is spin-dependent. It is shown that the interplay between the spin-dependent quasiparticle 
distribution and the triplet superconducting correlations induced by the proximity effect between 
the superconducting leads and ferromagnetic elements of the interlayer, leads to the appearence of an 
additional contribution to the Josephson current. This additional contribution jt can be extracted 
from the full Josephson current in experiment. The features of the additional supercurrent jt, which 
are of main physical interest are the following: (i) We propose the experimental setup, where the 
contributions given by the short-range (SRTC) and long-range (LRTC) components of triplet super- 
conducting correlations in the interlayer can be measured separately. It can be realized on the basis 
of S/N/F/N/S junction, where the interlayer is composed of two normal metal regions with a spiral 
ferromagnet layer sandwiched between them. For the case of tunnel junctions the measurement of 
jt in such a system can provide direct information about the energy-resolved anomalous Green's 
function components describing SRTC and LRTC. (ii) In some cases the exchange field-suppressed 
supercurrent can be not only recovered but also enhanced with respect to its value for non-magnetic 
junction with the same interface resistances by the presence of spin-dependent quasiparticle distri- 
bution. This effect is demonstrated for S/N/S junction with magnetic S/N interfaces. In addition, 
it is also found that under the considered conditions the dependence of the Josephson current on 
temperature can be nontrivial: at first the current rises upon temperature increasing and only after 
that starts to decline. 

PACS numbers: 74.45.-l-c, 74.50.-|-r 



I. INTRODUCTION 

Interplay between superconductivity and ferromag- 
netism in layered mesoscopic structures offers an arena 
of interesting physics to explore. By now it is already 
well known that so-called odd-frequency triplet pair- 
ing correlations are generated in hybrid superconduc- 
tor/ferromagnet (S/F) structures'^. The essence of this 
pairing state is the following. The wave function of a 
Cooper pair {ipcri ii'i,ti)ip(72 ('"2, ^2)) must be an odd func- 
tion with respect to permutations of the two electrons. 
Consequently, in the momentum representation the wave 
function of a triplet Cooper pair has to be an odd function 
of the orbital momentum for equal times ti = t2, that 
is, the orbital angular momentum L is an odd number. 
Thus, the triplet superconducting condensate is sensitive 
to the presence of impurities, because only the s-wave 
(L — 0) singlet condensate is not sensitive to the scat- 
tering by nonmagnetic impurities (Anderson theorem). 
S /F hybrid structures are usually composed of rather im- 
pure materials. Therefore, according to the Pauli princi- 
ple equal time triplet correlations should be suppressed 
there. However, another possibility for the triplet pairing 
exists. In the Matsubara representation the wave func- 
tion of a triplet pair can be an odd function of the Mat- 
subara frequency and an even function of the momentum. 
Then the sum over all frequencies is zero and therefore 
the Pauli principle for the equal-time wave function is 
not violated. These are the odd-frequency triplet pairing 
correlations, which are realized in S/F structures. 



If there is no a source of spin-flip processes in the 
considered structure (that is the magnetizations of all 
the magnetic elements, which are present in the sys- 
tem, are aligned with the only one axis) then the Cooper 
pairs penetrating into the nonsuperconducting part of the 
structure consist of electrons with opposite spins. Their 
wave function is the sum of a singlet component and a 
triplet component with zero total spin projection Sz = Q 
on the quantization axis. The resulting state has com- 
mon origin with the famous LOFF-state^ and can be 
reffered to as its mesoscopic analogue. This mesoscopic 
LOFF-state was predicted theoretically^i^ and observed 
experimentally2rii. In this state Cooper pair acquires the 
total momentum 2Q or —2Q inside the ferromagnet as a 
response to the energy difference between the two spin di- 
rections. Here Q cx h/vp, where h is an exchange energy 
and vp is the Fermi velocity. Combination of the two pos- 
sibilities results in the spatial oscillations of the conden- 
sate wave function "^{x) in the ferromagnet along the di- 
rection normal to the SF interfacei^. '^aix) (x cos(2(5a;) 
for the singlet Cooper pair and ^'t(x) cx sin(2Qx) for the 
triplet Cooper pair. The same picture is also valid in the 
diffusive limit. However, there is an extra decay of the 
condensate wave function due to scattering in this case. 
In the regime h ^ |A|, where A is a superconducting 
order parameter in the leads, the decay length is equal 
to the magnetic coherence length £^p = D jh, while 
the oscillation period is given by l-K^p. Here D is the 
diffusion constant in the ferromagnet, h — \ through- 
out the paper. Due to the fact that the decay length 
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is rather short (much less than the superconducting 
coherence length = y^D/ A) the sum of and 
(corresponding to Sz = 0) can be considered as a 
short-range component (SRC) of the pairing correlations 
induced by the proximity effect in the ferromagnet. 

The situation changes if the magnetization orienta- 
tion is not fixed. The examples are domain walls, spi- 
ral ferromagnets, spin-active interfaces, etc. In such a 
system not only the singlet and triplet Sz — com- 
ponents exist, but also the odd-frequency triplet com- 
ponent with 5z = ±1 arises in the nonsuperconducting 
region. The latter component penetrates the ferromag- 
net over a large distance, which can be of the order of 
= \/D/T in some cases. So, this triplet compo- 
nent can be considered as the long-range triplet com- 
ponent (LRTC). Various superconducting hybrid struc- 
tures, where LRTC can arise, were considered in the lit- 
erature (See Refs. 0, [H, [13 and references therein). In 
addition, the creation of LRTC was theoretically pre- 
dicted in structures containing domain walls'--^'-''^^, spin- 
active interfacesiiii^, spiral ferromagnets^-^-'^ and multi- 
layered SFS systems^i2^. There are several experimen- 
tal works, where the long-range Josephson effect^"— and 
the conductance of a spiral ferromagnet attached to two 
superconductors21 were measured. These results give 
quite convincing evidence of LRTC existence. 

Practically all the discussed above papers are devoted 
to investigation of an odd-frequency triplet component 
under the condition that the energy distribution of quasi- 
particles is equilibrium and spin-independent. However, 
as it was shown recently^, the creation of spin-dependent 
quasiparticle distribution in the interlayer of SFS junc- 
tion leads to appearence of the additional contribution 
to the Josephson current through the junction. This 
additional supercurrent flows via vector part AT,- 1 of 
supercurrent-carrying density of states, which does not 
contribute to the Josephson current in a junction with 
s-wave superconductor leads if the quasiparticle distri- 
bution in the interlayer is spin-independent. Below we 
briefly describe how this effect arises. 

The energy spectrum of the superconducting correla- 
tions is expressed in a so-called supercurrent- carrying 
density of states (SCDOS)22r— . This quantity represents 
the density of states weighted by a factor proportional to 
the current that each state carries in a certain direction. 
Under equilibrium conditions the supercurrent can be ex- 
pressed via the SCDOS as^i 

j oc J deNj{e)ta.nhe/2T , (1) 

where e stands for the quasiparticle energy, tanhe/2r ~ 
ip{s) is the equilibrium distribution function and Nj{e) is 
SCDOS. In the presence of spin effects SCDOS becomes 
a matrix 2 x 2 in spin space and can be represented as 
Nj — Nj^s + Nj^txr, where ai are Pauli matrices in spin 
space. Scalar in spin space part of SCDOS Nj^s is referred 
to as the singlet part of SCDOS in the paper and vector 
part Nj,t is referred to as the triplet part. Nj,t is directly 



proportional to the triplet part of the condensate wave 
function. It is well known that the spin supercurrent 
cannot flow through the singlet superconducting leads. 
Therefore, Nj^t does not contribute to the supercurrent 
in equilibrium. Having in mind that the triplet part of 
SCDOS is even function of quasiparticle energy, one can 
directly see that this is indeed the case. Otherwise, if 
the distribution function becomes spin-dependent, that 
is (p{e) = <(5o(£) + ¥'(£)"■! the supercurrent carried by 
the SCDOS triplet component Nj^t in the ferromagnet 
is non-zero because the scalar product Nj^t{£)v{£) con- 
tributes to the spinless supercurrent in this case^^. 

As it is obvious from what discussed above, the spin in- 
dependent nonequilibrium quasiparticle distribution does 
not result in an additional contribution to the supercur- 
rent flowing via Nj^t- However, it is worth noting here 
that the effect of the spin independent nonequilibrium 
distribution function has been considered as well^^i^. 
It was shown that in the limit of small exchange fields 

<C I A I the combined effect of the exhange field and the 
nonequilibrium distribution function is also nontrivial. 
For instance, part of the field-suppressed supercurrent 
can be recovered by adjusting a voltage between addi- 
tional electrodes, which controls the distribution func- 
tion. 

In the present paper we continue investigation of 
the interplay between the triplet correlations and spin- 
dependent quasiparticle distribution. As it was explained 
above, the simultaneous presence of the triplet correla- 
tions and spin-dependent quasiparticle distribution in the 
interlayer results in appearence of the additional con- 
tribution to the supercurrent flowing via Nj^t- In the 
present paper we concentrate on two features of this ad- 
ditional supercurrent, which are of main physical interest 
and propose appropriate mesoscopic systems, where they 
can be observed: 

(i) The additional supercurrent allows for direct mea- 
surement of the energy-resolved odd-frequency triplet 
anomalous Green's function in the interlayer. The point 
is that for junctions with low-transparency interfaces be- 
tween the superconductor and the interlayer region Nj^t 
is directly proportional to the triplet part of anomalous 
Green's function in the interlayer. By measuring the 
"non-local" conductance (that is, the derivative of the 
critical current with respect to voltage V, which is ap- 
plied to the additional electrodes attached to the inter- 
layer region and controls the value of spin injection into 
the interlayer), one can experimentally obtain the value 
of triplet part of the anomalous Green's function in the 
interlayer as function of energy. As it was discussed in 
the introduction, the triplet correlations induced by the 
proximity effect in S/F structures are odd in Matsub- 
ara frequency, that is the corresponding two-particle con- 
densate wave function taken at coinciding times is zero. 
Therefore, the direct measurement of the energy-resolved 
anomalous Green's function is of great interest. 

Here we propose an experimental setup, which allows 
for extracting from the current SRTC and LRTC con- 
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tributions and their separate observation. By measuring 
the "non-local" conductance one can separately obtain 
the values of LRTC and SRTC triplet parts of the anoma- 
lous Green's function in the interlayer as functions of en- 
ergy. It is based on a multilayered S/N/F/N/S junction, 
where a layer made of a weak ferromagnetic alloy having 
exchange field A<</i<<£i?is sandwiched between 
two normal metal layers. The direction of the F layer 
magnetization is assumed to be nonuniform in order to 
have a possibility of LRTC investigation. The leads are 
made of dirty s-wave superconductors. 

While all the experiments described in the introduc- 
tion give unambigous signatures of the fact that the odd- 
frequency triplet correlations do exist in hybrid SF sys- 
tems, they do not allow for direct investigation of the 
triplet anomalous Green function in dependence on en- 
ergy. For example, Josephson current in equilibrium is 
only carried by the scalar part of SCDOS Nj^s- Surely, 
it is modified by the presence of the triplet component 
(and, in particular, manifests weakly decaying behavior if 
LRTC is present in the system). However, Nj^g is not di- 
rectly proportional to the triplet anomalous Green func- 
tion, but can contain it only in a nonlinear way. The 
other measurable quantity in equilibrium is the local den- 
sity of states (LDOS), where the odd- frequency triplet 
component manifests itself as a zero-energy peak. This 
effect has been studied as for SF bilayers so as for SN 
bilayers with magnetic interfaces'^"—. However, LDOS 
is also not directly proportional to the triplet anomalous 
Green's function. The oscillating behavior of the crit- 
ical temperature as a function of an SF bilayer width 
(See, for example, Ref. TT and references therein) is also 
an excellent fingerprint of the triplet correlations (one- 
dimensional LOFF state) presence. However, the or- 
der parameter in the singlet superconductor S is related 
only to the singlet part of the anomalous Green func- 
tion, which is modified by the presence of triplet correla- 
tions, but does not allow for their direct observation. On 
the other hand, in case if the quasiparticle distribution 
is spin-dependent, quantities, which are directly propor- 
tional to the triplet anomalous Green function start to 
contribute to experimentally observable things. Joseph- 
son current under the condition of spin-dependent quasi- 
particle distribution in the interlayer is one of them. 

(ii) It is well-known that fcrromagnetism and singlet 
superconductivity are antagonistic to each other. In over- 
whelming majority of situations it results in the suppres- 
sion of the Josephson current through the system with 
ferromagnetic elements with respect to the system with 
the same interface resistances but without ferromagnetic 
elements. This is also valid even if LRTC is formed in 
the system. In the present paper we show that in some 
cases the exchange field-suppressed supercurrent can be 
not only recovered but also enhanced with respect to its 
value for non-magnetic junction with the same interface 
resistances by the presence of spin-dependent quasiparti- 
cle distribution. That is, roughly speaking, in some cases 
the spin-dependent quasiparticle distribution can over- 



compensate the suppression of proximity-induced super- 
conducting correlations by fcrromagnetism. We demon- 
strate that such an effect can be observed in S/N/S junc- 
tion with magnetic interfaces. 

The paper is organized as follows. In Sec. [Tl] the con- 
sidered model systems are described and the theoretical 
framework to be used for obtaining our results is estab- 
lished, in Sec. IIIII we present the results of the Joseph- 
son current calculation for a multilayered S/NFN/S sys- 
tem under spin-dependent quasiparticle distribution and 
demonstrate how to obtain from these data information 
about the structure of the odd-frequency triplet correla- 
tions. Sec.|TV]is devoted to consideration of SNS junction 
with magnetic SN interfaces under similar conditions for 
the quasiparticle distribution in the interlayer. We sum- 
marize our finding in Sec. |Vl In Appendix |^ we rep- 
resent the results for anomalous Green's function in the 
interlayer and all the parts of the Josephson current for 
S/NFN/S junction, calculated in the framework of a par- 
ticular microscopic model of N/F/N layer. Appendix [B] 
is devoted to a particular microscopic model of magnetic 
S/N interface, which we assume to be more appropri- 
ate for the investigation of current enhancement in SNS 
junction. 

II. MODEL AND GENERAL SCHEME OF 
CALCULATIONS 
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FIG. 1. S/N/F/N/S junction under consideration with the 
additional electrodes, which are proposed to be used for cre- 
ation of a spin-dependent quasiparticle distribution in the in- 
terlayer. 

The first system we consider is a multilayer S/NFN/S 
Josephson junction shown schematically in Fig. [1] It 
consists of two s-wave superconductors (S) and an in- 
terlayer composed of two normal layers iV' and iV with 
a ferromagnetic layer F, sandwiched between them. The 
X-axis is directed along the normal to the junction and 
the y- and z-axes are in the junction plane. The co- 
ordinates of FN interfaces are x = ^dp/^, while SN 
interfaces are located at x = ^{dp + dN)/2. That is. 
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the full length of the F layer is dp, while the length 
of each N layer is 6,^/2. The middle F layer is sup- 
posed to have the exchange field h satisfying the con- 
dition Is. ^ h £p. The exchange field of F layer is 
assumed to be non-homogeneous, what allows for the ex- 
istence as triplet pairs with opposite spins (SRTC) so as 
triplet pairs with parallel spins (LRTC) in the interlayer. 
h = /i(0, sin0(a;), cos0(x)), that is the magnetization 
vector rotates in the F layer (within the junction plane) . 
For simplicity we suppose that the rotation angle has a 
simple x-dependence 



e(x)=e'x, -dpl2<x<dFl2 



(2) 



where 0' does not depend on coordinates. The additional 
electrodes are supposed to be attached to the N layers 
in order to make it possible to create a spin-dependent 
quasiparticle distribution in the interlayer. 

We use the formalism of quasiclassical Green-Keldysh 
functions^ and assume that the superconductors and 
all the internal layers are in the diffusive regime. The 
fundamental quantity for diffusive transport is the mo- 
mentum average of the quasiclassical Green's function 
g{r,e,t) = {g{pf,r,e,t))p^. It is a 8 x 8 matrix form 
in the product space of Keldysh, particle-hole and spin 
variables. In the absence of an explicit dependence on 
time variable the Green's function g{r,e) in the inter- 
layer obeys the Usadel equation 



D 



— V(5V5) + [ersaopo -h,g] =0 



(3) 



where Ti, ai and pi are Pauli matrices in particle-hole, 
spin and Keldysh spaces, respectively, tq, ctq and po 
stand for the corresponding identity matrices. For sim- 
plicity, the diffusion constant D is supposed to be iden- 
tical in all three internal layers. The matrix structure of 
the exchange field is as follows 



h = h(Tpa{l + T3)/2 + h(T*poil - T3)/2 



(4) 



The exchange field h rotates in the F layer according to 
the described above model. In the N layers h — 0. 

The Usadel equation ([3]) should be supplied with the 
normalization condition g^ = —■n^TQU^pQ and is subject 
to Kupriyanov-Lukichev boundary conditions'^ at S/N 
and N/F interfaces. The barrier conductances of the left 
and right S/N interfaces are assumed to be identical for 
simplicity and are denoted as Gt- Then the boundary 
conditions at S/N interfaces take the form 



(jNdxgN 



Gt 
— 

2a N 



[9N,gs] 



(5) 



Here ^at is the solution of the Usadel equation ^ at the 
left {x = -{dN + dF)/2) or right {x = {dN + dF)/2) S/N 
interface, a — -1-1 (—1) at the left (right) interface, cttv 
is the conductivity of N layers and ap is the conductiv- 
ity of F layer (defined for later use), gs stands for the 
Green's functions at the superconducting leads. Due to 



the fact that we are mostly interested in the case of low- 
transparent S/N interfaces below, we can safely neglect 
the suppression of the superconducting order parameter 
in the S leads near the interface and take the Green's 
functions at the superconducting side of the boundaries 
to be equilibrium and equal to their bulk values. In this 
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where k = +1(— 1) for the retarded (advanced) Green's 
function, x stands for the order parameter phase differ- 
ence between the superconducting leads and (5 is a posi- 
tive infinitesimal. 

The second model system, which we consider in or- 
der to study the enhancement of field-suppressed su- 
percurrent under the spin-dependent distribution, is an 
S/N/S junction with magnetic S/N interfaces. The full 
length of the normal region is djv, the x-axis is normal 
to the junction plane and the interfaces are located at 
X = =pdjv/2. As in the previous case, additional elec- 
trodes are attached to the interlayer region for creation of 
a spin-dependent quasiparticle distribution in the inter- 
layer. The Green's function in the N layer obeys Eq. ([3]) 
provided that h = 0. However, the boundary conditions 
contain additional terms with respect to Eq. JS]) because 
the transmission properties of spin-up and spin-down 
electrons into a ferromagnetic metal or a ferromagnetic 
insulator are different, which gives rise to spin-dependent 
conductivities (spin-filtering) and spin-dependent phase 
shifts (spin-mixing) at the interface. The generalized 
boundary conditions for the diffusive limit can be written 
in the formiSiii 

mdxgN = -a- — [gN,gs] - 
ZaN 

-a^^ [ffAT, {rha,gs}] + [mQ,^^] , (9) 
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where g^ is the Green's function value at the normal 
side of the appropriate S/N interface (at x = =FdAr/2). 
As above, gs stands for the Green's function in the 
superconducting lead and is expressed by Eqs. ([S])-®. 
rha = ma(Tpa{l + t^)/2 + maa*po{l - t^)/2, where 
rria is the unit vector aligned with the direction of the 
left [a = +1) or right [a — —1) SN interface magneti- 
zation. {...} means anticommutator. The second term 
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accounts for the different conductances of different spin 
directions and Gmr ~ Gt,^ — Gt,i- The third term 
~ G^ gives rise to spin-dependent phase shifts of quasi- 
particles being reflected at the interface. It is worth ot 
note here that boundary conditions (|9]) are only valid for 
small (with respect to unity) values of transparency and 
spin-dependent phase shift in one transmission channel^. 
However, for the case of plane diffusive junctions we can 
safely consider — G^£_s/<^n > 1 due to a large number 
of channels. G^ has been calculated for some particular 
microscopic models of the interface^Sii^ and can be large 
enough even if the conductance Gt — >■ 0. In Appendix 
|B]we calculate G^, Gt and Gmr for S/N interface com- 
posed of an insulating barrier and a thin layer of a weak 
ferromagnetic alloy. We suppose this microscopic model 
to be the most appropriate to the considered problem. 

In what follows we assume that the S /N interfaces are 
low-transparent for the both considered systems, that is 
Gt = Gt^s/o-n < 1 and Gmr = Gmr^s/o^n < 1- 
In order to calculate the Josephson current through the 
junction in the leading order of the interface transparency 
Gt it is enough to obtain the retarded and advanced 
Green's functions in the leading order of the trans- 
parency. If one makes use of the following definitions for 
the Green's function elements in the particle-hole space 
(all the matrices denoted by ' are 2x2 matrices in spin 
space throughout the paper) 



/ 'R.A fR,A 



(10) 



then one can obtain from the Usadel equation ([3]) and 
the appropriate boundary conditions Eq. ([5]) or ([9]) that 
the diagonal in particle-hole space elements of g^'"^ are 
zero-order in Gt quantities and take the following form 
in the interlayer 



linearized Usadel equations for the anomalous Green's 
function ^f^'^, //^'"*| can be written as follows 

cR-A 



R,A 



1 rR.A 



2£/«-- ~ 2/i/f - iKDdtn^^- = (13) 
According to the general symmetry relation^ f^'^i^) — 
f^'^*{—e) the singlet and triplet parts of /^'"^ can be 
expressed via the corresponding parts of /^'^ as follows 



/,(£) = -/;(-£) 



(14) 



The linearized Usadel equations dTS]) should be sup- 
plemented by the appropriate boundary condition, which 
are to be obtained by linearization of Eq. ©or © and 
at the S/N interfaces take the form 
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where f^'f and f^'f are the singlet and triplet part 
values of the anomalous Green's function at the normal 
side of the S/N interface. G^ ^ only if the S/N interface 
is spin-active. It is worth to note here that in this linear 
in Gt and Gmr approximation the term proportional to 
Gmr does not enter the boundary conditions. 

Eqs. and ([15]) allow for the calculation of the re- 
tarded and advanced Green's functions in the leading in 
transparency approximation. However, it is not enough 
for obtaining of the electric current through the junction, 
which should be calculated via Keldysh part of the qua- 
siclassical Green's function. For the plane diffusive junc- 
tion the corresponding expression for the current density 
reads as follows 



^R,A ^ 



-\-oo 



(11) 



The off-diagonal in particle-hole space elements of the 
Green's function are of the first order in Gt and should be 
obtained from the linearized Usadel equations, which are 
to be derived from Eq. ([3]). It is convinient to represent 
the off-diagonal elements in the following form 



rR,A ■ 
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where f^'^ and f^'^ denote the singlet and triplet parts 
of the anomalous Green's function, respectively. For the 
case we consider (the magnetization vectors of all the 
ferromagnetic layers and spin-active interfaces, which are 
present in the system, are in the junction plane) the out- 
of-plane a;-component of the triplet part is absent and the 
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— {9{x,e)d:,g{x,e))^ 



(16) 

where e is the electron charge. The expression is writ- 
ten for the normal layer, but it is also valid for the 
ferromagnetic region with the substitution ap for (Tat. 
{g{x,e)dxg{x,e))^ is 4 x 4 Keldysh part of the corre- 
sponding combination of full Green's function. It is con- 
venient to calculate the current at the S/N interfaces. 
Then the required combination of the Green's functions 
can be easily found from Keldysh part of boundary con- 
ditions ([5|) or ([5]) . In addition, we express Keldysh part of 
the full Green's function via the retarded and advanced 
components and the distribution function: g^ = g^(f — 
<pg^. Here argument {x,e) of all the functions is omit- 
ted for brevity. The distribution function is diagonal in 
particle-hole space: cp = ip{To + T3)/2 + a2'fcr2{To — T3)/2. 
Then to the leading (second) order in transparency 
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2i7r (cosh Q§ + cosh 65) ( tanh 



2T 



2i7r (cosh 9f + cosh 65) ma^p 



(o)+(i) 



, (17) 



where f^'f and /j^'^ are taken at the normal side of the 
appropriate S/N boundary, (po and <p represent the scalar 
and vector parts of the distribution function ip = ipo+ifcr, 
which is also taken at the normal side of the appropriate 
S/N boundary. The superscripts and of 

the distribution functions mean that the corresponding 
quantity is calculated up to the zero and the first orders of 
magnitude in the interface conductance Gt, respectively. 

In order to calculate the current through the junction 
one should substitute Eq. p7)) into Eq. ([T6| . The result- 
ing expression for the current can be further simplified 
by taking into account the general symmetry relations 
between the Green's function elements^ expressed by 
Eq. (fH|) and the ones given below 

ffis) = m-e) 

<^o(e) = -<y5o(-£) 
cp(e) = ^{-e) . (18) 

Then the expression for the current density takes the 
form 
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the case of spin-independent distribution function. They 
do contribute, as it was demonstrated by a number of 
experiments discussed in the introduction. The point is 
that for the considered case of the tunnel junction fs in 
general contains long-range contributions resulting from 
the LRTC (if they are present in the system). It is worth 
to emphasize that all the aforesaid only concerns the cal- 
culation of the current at the interface. If one would cal- 
culate the current at an arbitrary point of the interlayer, 
the corresponding expression would contain ft quadrat- 
ically. Surely, the current by itself does not depend on 
x-coordinate, as it is required by the current conserva- 
tion. 



The distribution function entering current (fT9)) 

should be calculated by making use of the kinetic equa- 
tion, which is obtained from the Keldysh part of Usadel 
equation ([3]). As we only need the distribution function 
up to the first order in the interface conductance, all the 
terms accounting for the proximity effect (which are of 
the second order in Gt) drop out and the kinetic equa- 
tion takes especially simple form (we do not take into 
account inelastic relaxation in the interlayer) 



cR iax/2 



tanh— Re [sinhef] 



Re 
Re 



R „iaxl1 



N,. 



fR iax/2 



p>^^hxn [sinhef] + 
(^(•'^Im [sinhef] - 
.coshef[4"^+(^'(e) + 4"^+'^^(-£)]/2) + 
uGmr^ cosh ef m„ [vp(o)+(i) (e) + (-e)] } . (19) 

The additional contribution to the current, which is ab- 
sent for a spin-independent distribution function, is given 
by the third term. As it was mentioned in the introduc- 
tion, this term (connected to the triplet part of SCDOS) 
is directly proportional to the triplet anomalous Green's 
function at the interface. The fifth term also results from 
vector part of the distribution function, but under the 
considered conditions it does not contribute to the cur- 
rent, as it is shown below. It is worth to note here that, 
as it is seen from Eq. (|T9)) . the singlet part of SCDOS 
is only expressed via the singlet part of the anomalous 
Green's function. However, it does not mean that the 
triplet correlations do not contribute to the current for 



— [h{x)(r,pi\ = 



(20) 



where the exchange field h{x) is determined above for 
the ferromagnetic layer and vanishes for all the normal 
regions. 



The kinetic equation should be supplemented by the 
boundary conditions at the S/N interfaces and the in- 
terfaces with additional electrodes, attached to the nor- 
mal regions of the interlayer in order to create a spin- 
dependent quasiparticle distribution. While the bound- 
ary conditions at the interfaces with additional electrodes 
are discussed below for a particular considered system, 
the boundary conditions at the S/N interfaces are ob- 
tained from the Keldysh part of Eqs. (jH]) or ([S]) and to 



7 



the first order in the interface conductance take the form 



"^(coshef 



2a 



N 



niacr, (fi 



fsR is rapidly decaying in the interlayer, while Jlr 
is slowly decaying. Let us consider fsR at the left S/N 
interface (the left interface is chosen just for definiteness) . 
It can be rewritten in the form 



2(7 AT 

aGMR 

UN 



coshe4) (tanh 



e 
e 



(0) 



(cosh 85 + cosh 85) tanh ■^nriaCr 



coshOgniaCrip''^^ — cosh Qgip^^^m aO- 
aiG, 



2aN 



niacr, If 



(1) 



SR 



SR" 



/; 



(23) 



where /^^ is generated by the proximity effect at the 
left S/N interface itself and /^^ comes from the right 
S/N interface. It can be shown that for thick enough 
ferromagnetic layer dp/S^p » 1 fgn/fsi^ is proportional 
the small factor e '^''Z^'' . On the contrary, if Jlr is 
represented as 



For the case of multilayered N/F/N interlayer Eq. (|20|) 
should be also supplemented by boundary conditions 
at the N/F interfaces, which are to be obtained from 
Keldysh part of Kupriyanov-Lukichev boundary condi- 
tions ^ and are given in Appendix |X] for a particular 
microscopic model. 



III. S/NFN/S JUNCTION 

Now we consider the particular systems. This section is 
devoted to S/NFN/S Josephson junction. The model as- 
sumed for the exchange field of the F layer is already de- 
scribed above. The anomalous Green's function is found 
up to the first order in S/N conductance Gt according to 
Eqs. ([12]), and boundary conditions at F/N interface, 
which should be easily obtained from Eq. ([5]) for a given 
conductance of this interface. We assume that the mag- 
netization of the F layer rotates slowly, that is Q'£,p <C 1, 
while 8'^s ^ 1 or even larger than unity. This assump- 
tion seems to be quite reasonable^. Therefore, upon cal- 
culating the anomalous Green's functions we disregard 
all the terms proportional to Q' ^JD/h and higher pow- 
ers of this parameter, while keeping the terms, where 8' 
enters in the dimensionless combination Q' ^/D/\e\. 

To this accuracy the triplet part of the anomalous 
Green's function can be represented as 

ft' = (0, /„ /.) , 
fy = smQ{x)fsR{x) - cosQ{x) fLR{x) , 

= cos eix)fsR{x) + sin e{x)fLRix) , (22) 

where the z-axis is aligned with the direction of the ex- 
change field in the middle of the F layer (at a; = 0). fsR 
Ulr) is formed by the Cooper pairs composed of the elec- 
trons with opposite (parallel) spins. We are interested 
in the values of the triplet component at the S /N inter- 
faces, where sin8(a;) = — asin8 = — asin [8'(ii?/2]. The 
particular expressions for Jsr, and Jlr and singlet com- 
ponent fs depend strongly on the conductance of F/N 
interface and in general are quite cumbersome. In order 
to give an idea of their characteristic behavior we have 
calculated them for the most simple model of absolutely 
transparent F/N interfaces. The corresponding expres- 
sions are given in Appendix VK\ 



Jlr. — f\ 



LR^ 



(24) 



fbR does not contain the small factor e~'^''/^'^ in the 
leading approximation and, therefore, Jlr describes the 
LRTC. As it is explicitly demonstrated in Appendix \K[ 
the characteristic decay length of Jlr in the F layer is 

|At|-\ where A* = ^^8'^ - 2i{e + iS)/D. It is much 

larger than £^p for the considered case <C 8' ^. 

To the considered accuracy the singlet component of 
the anomalous Green's function also decays at the dis- 
tance ~ in the F layer, just as SRTC Jsr does, be- 
cause it is also composed of the electron pairs with an- 
tiparallel spin directions. Indeed, if at the left bound- 
ary is also represented as 



(25) 



then /J cx e'"^"/^'' in the limit dp/^p > 1. 

Therefore, the main contribution to the Josephson cur- 
rent Eq. (fT9|) is given by the LRTC component /lr of the 
anomalous Green's function. However, this contribution 
is nonzero only for the case of spin-dependent quasipar- 
ticle distribution. In the standard case of thermal spin- 
independent quasiparticle distribution the current is de- 
termined by the singlet component /g. Consequently, it 
only contains the term proportional to the small factor 
Q-dp/iF ^ p^^ fjj-g^ glance, it contradicts to the well-known 
fact that the equilibrium Josephson current contains the 
contribution generated by the LRTC, if it is present in 
the systemi*^. In fact, if one calculates the current at 
the S/N boundary, then should be modified by pres- 
ence of LRTC and should contain a slowly decaying term, 
which provides the appropriate contribution. It is indeed 
the case for the system we consider. However, the cor- 
responding term is proportional to (8'^f)^ and is dis- 
regarded in our calculation. Surely, it should be taken 
into account upon calculating the Josephson current for 
the case of spin-independent quasiparticle distribution, 
because in spite of the small factor {Q'£,pY it can result 
in large enough current contribution due to the absence 
of the suppression factor e"''^/^^. At the same time we 
can safely disregard this term, because for the consid- 
ered case of spin-dependent quasiparticle distribution the 
main contribution to the Josephson current is given by 
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/lr term, which contains neither (Q'^p)"^ nor the ferro- 
magnetic suppression factor e"'^''/^" . 

In order to generate a spin- dependent quasiparticle dis- 
tribution in the interlayer, additional electrodes are at- 
tached to the normal regions of the interlayer. While in 
the paper we propose some particular way of such a dis- 
tribution creation, it is not important how particularly it 
is obtained. The main point is to have a vector part 
of the distribution function in the interlayer, generated 
anyway. For example, it can be created by a spin injec- 
tion into the interlayer. If this is the case, the discussed 
below results qualitatively survive. 

In the present paper we assume that each of the nor- 
mal regions of the interlayer is attached to two additional 
normal electrodes N^^ and (See Fig. In their 

turn, the electrodes Nl and Nl have insertions Pi and 
Pr made of a strongly ferromagnetic material. Let the 
F''^''^ is applied between the elec- 



voltage V/^''^ - Vl^' > = 
trodes Ni^""^ and nI^'' 



Here V^ 



'^'^^ and v/'''"'' are the 

l(r) 



electric potentials of the outer regions of the N\; ' and 
electrodes with respect to the potential of the su- 



n: 



:(r) 



perconducting leads. It is worth to note here that the 
superconductor is assumed to be closed to a loop and 
the voltage between the superconducting leads is absent. 
The conductances of the 7V'('')/iV^'''^ and N^^^'^nI^"'^ in- 
terfaces are denoted by g^^^^ and g['"^\ respectively. 

Further, for definiteness we consider the left normal re- 
gion of the interlayer with the corresponding additional 
electrodes. We choose the quantization axis z; along the 
magnetization of the left ferromagnetic insertion Pi and 
the definitions i?p,-|-, Rpii stand for the Pi region resis- 
tivities for spin-up and spin-down electrons. Then under 
the conditions that (i) the iV' layer resistance Rn and 
the resistance of region inclosed between iV' and Pi 
can be disregarded as compared to l/gt and Rpii and 
(ii) 1/Rp,i <^ g\-€. 1/Rpif one can believe that the volt- 
age drops mainly at the Pi region for spin-down electrons 
and at the N'-fN^ interface for spin- up electrons. Also, 
the dissipative current flowing through Nl/N^ /Nl sys- 
tem is small and can be disregarded. Consequently, it 
is obtained that the electric potentials for spin-up and 
spin-down electrons in the iV^ region inclosed between Pi 
and are different and practically constant over this 
region. While the spin-up electrons are at the electric 
potential VJ,' in this region, the potential for spin-down 
electrons is approximately V} . 

In order to simplify the calculations we assume that 
= V'' = Vb-Vt = 2V. The left and the right addi- 
tional electrodes only differ by the direction of the mag- 
netization of the Pi and Pr insertions. For later use we 
define the unit vectors aligned with the Pi and Pr magne- 
tizations as Ml and iV/r, respectively. In order to satisfy 
the electroneutrality condition the electric potential of 
the superconducting leads should be equal to {Vt + Vh)/2. 
Then the electric potentials for spin-up and spin-down 
electrons in the Nl region inclosed between Pi and N^ 



counted from the level of the superconducting leads are 
Ft = {Vb - Vt)/2 = V and ^ = {Vt - Vb)/2 = -V . 
Due to the fact that one can disregard the voltage drop 
inside this region, the distribution functions for spin-up 
and spin-down electrons in this region are close to the 
equilibrium form (with different electrochemical poten- 
tials) . For the general case (if the quantization axis does 
not aligned with the Pi magnetization) the distribution 
function becomes a matrix in spin space and takes the 
form 



= <PoCTo + (fitMia , 



tanh ■ 



eV 



tanh 



2T 

e-eV 
2T 



tanh 



e + eV 



tanh 



2T 

e + eV 
2T 



(26) 



The same form of the distribution function is valid for 
the NJ^ part enclosed between Pr and N"^ with the sub- 
stitution Mr for Ml. 

Now we can obtain the distribution function in the 
iV' and TV regions of the interlayer, which enters cur- 
rent (fT9|) . For the considered case gt ^ 1 the dissipative 
current flowing through 7V^/7V'/iV/ junction is negligible 
and, therefore, the y-dependence of the distribution func- 
tion in the iV'*^''-' region can be disregarded. Then under 
the condition cjp aj^ the distribution functions 
m the N^ and TV regions calculated up to the zero order 
in the S/N conductance Gt are spatially constant and 
equal to ipi and Lpr^ respectively. Indeed, the equation 
for the distribution function ((20|) aX h — Q and bound- 
ary conditions at S/N interfaces (pij) (corresponding to 
= 0) are satisfied by this solution. Boundary condi- 
tions at F/N interfaces (|A3|) are satisfied approximately 
due to the smallness of the distribution function gradient 
under the condition ap <Si crjv- If one goes beyond the 
approximation apjaj^ <C 1, the distribution function in 
the TV' and W regions acquires gradient terms propor- 
tional to the parameter (Jp/(Jn- If the F/N interface is 
less transparent than it is considered in Appendix |3 the 
distribution function gradient in the N layer even smaller 
and the condition tii? <C (Jn is not so necessary. 

Although the distribution function in the middle F 
layer does not enter the current expression ([T9|) . it is 
interesting to discuss here how it behaves. For simplicity 
we consider the limiting case Q' — >■ 0, when the exchange 
field in the ferromagnet is practically constant and the 
qualitative physical picture is more clear. According to 
Eq. ((20|) and boundary conditions (IA3P the scalar part of 
the distribution function (^o is constant over the F layer 
and coincide with its value in the iV' and N"^ regions. The 
vector component parallel to the exchange field of the 
ferromagnet is a linear function of x-coordinate, which 
matches the constant values iptMi(^r)hi(r}/h at the F/N 
interfaces. Here hi r = h{x = ^dp/2) are the exchange 
field values at the left and right N/F interfaces. The 
vector component perpendicular to the exchange field of 
the ferromagnet decays from the F/N interfaces into the 
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ferromagnetic region at the characteristic length os- 
cillating simultaneously with a period , as it can be 
obtained from Eq. (|20|) and boundary conditions ()A3|) . 

Strictly speaking, the distribution function in the N 
layers only takes form (|26p if one assumes no spin re- 
laxation there. Spin relaxation processes reduce vector 
part ipt of distribution function (j26p . The reduction can 
be roughly estimated as ipf^ = ipt/{l + Tesc/Tsr)- Here 
ipf^ is the vector part of the distribution function in the 
presence of spin relaxation processes, while (pt is defined 
by Eq. (f26|) . Tesc = o-i^dy/D{gb + gt) is an effective time, 
which an electron spends in the N layer before escaping. 

is the characteristic spin relaxation time. So, spin re- 
laxation processes do not qualitatevely influence on the 
distribution fuction if Tesc/Tsr ^ 1- This condition seems 
to be not restrictive in real materials. For example, if 
one assume that the N layers are made of Al in normal 
state, where Xgr ~ \/ Dtst — 450//to'*- has been reported 
and the condition t^sc/tst — <^Ndy/{gb + gt)^^sr ^ 1 '^^'^ 
be valid in a wide range of the values of dimensionless 
parameter [gi, -\- gt)^s/<^N characterizing the joint con- 
ductance of TV^'^'Va^''''^ and nI'^^'^/N^'^^'I interfaces. 

Now we turn to the discussion of the Josephson cur- 
rent through the junction. It is expressed by Eq. (1191) . 
As for the considered case of nonmagnetic S /N interfaces 
Gmb. = 0, the last term in this formula is absent. Due 
to the fact that the scalar part ip'j^'' of the distribution 
function in the interlayer [Eq. (pS)) ] is an odd function 
of quasiparticle energy, the part of the current generated 



and (jlSp ] this contribution takes the form 



■ Lr 



hl.rM, 



[Mur X hur) ex 



by the term c>c cosh 9^ 



/2 also van- 



ishes. Further, in order to avoid the flowing of a quasi- 
particle current through the junction we assume that 
\eV\ < A and the temperature is low (T <C A). Under 
these conditions the linear in x-coordinate part of ip^^"^ 
(it is this term that provides the flowing of the quasipar- 
ticle current through the junction) is zero in each of the 
N regions of the interlayer, as it is dictated by boundary 
conditions (|2ip . Therefore, (^*^^-' is approximately con- 
stant in each of the N layers. We comment on the values 
of these constants below. 

The first two terms in Eq. p9l) represent the contribu- 
tion of the SCDOS singlet part, which takes place as for 
the case of spin-independent quasiparticle distribution, 
so as when this distribution is spin-dependent. We re- 
fer to this contribution as jg. The particular expressions 
for js can be easily found in the framework of a given 
microscopic model of NFN interlayer after substitution 
of the particular expressions for the singlet part of the 
anomalous Green's function and the scalar part of the 
distribution function [Eqs. ((261) and ([18])] into Eq. (fT9)) . 

The third term in Eq. (|19p contains the current flow- 
ing through the SCDOS triplet part and is nonzero only 
for the case of spin-dependent quasiparticle distribution. 
This contribution is the main result of the present sec- 
tion. If one substitutes the particular expressions for the 
triplet part of the anomalous Green's function [Eq. (P^ 
and the vector part of the distribution function [Eqs. (pS)) 



Jt = -JSR ^ ' + ajLB- — — ^ ' ' , (27) 

where is the unit vector along the x-direction. The 
currents jsB. and jLi?, are generated by the short-range 
and long-range triplet components of the anomalous 
Green's function, respectively. Consequently, if the F 
layer is thick, that is £,f ^ dp, the current jsR (as well 
as js) is small due to the factor e"''^/^^, while Jlr is 
not suppressed by this factor. The particular expression 
for js, jsR and j^n calculated in the framework of the 
considered here microscopic model are given in Appendix 

El 

It is seen from Eq. ([TT)) that the values of the current 
contribution jt at the left and right S/N interfaces can 
be different, that is, in general, jl'^ = jt ± ja- However, 
under the condition that the superconducting leads are 
closed into a loop the currents at the left and right S/N 
interfaces must be equal to each other. It appears that 
the distribution function in the N layers acquires addi- 
tional terms V'/ which are proportional to Gt- Under 
the condition V'' = y we obtain that (/sp' = (pi^^ . Then, 



according to Eq. (|T9l) this term results in the current con- 
tribution, which exactly compensates ja- Therefore, the 
Josephson current jt flowing through the junction can be 
simply calculated as jt = (j| -I- jt)/'^- 

It is obvious from Eq. (j26p that js is an even function 
of voltage V applied to the additional electrodes and jt 
is an odd function of this voltage. Therefore, it is easy 
to extract in experiment contributions js and jt from 
the full Josephson current: js{V) = ij{V) + j{~V))/2, 
while jt{V) = {j{V) — j{—V))/2. Further, it is seen from 
Eq. (|27p that by choosing the appropriate orientation of 
Pi and Pr magnetizations, one can, in principle, measure 
either jsR or jpn current contributions. For this reason 
it makes sense to discuss all the current contributions js , 
jsR. and jLR separately. In the tunnel limit all of them 
manifest sinusoidal dependence on the superconducting 
phase difference x, that is js.sR.LR = Jlsfl.Lfl^inX- 
Therefore, we only discuss the corresponding critical cur- 
rents , j%jf , jlji below. 

Fig. [2] represents these contributions, calculated in the 
framework of the microscopic model discussed in Ap- 
pendix El as a function of voltage V- First of all, it 
is worth to note that current components jsR^ and Jla, 
carried by the triplet part of SCDOS, are nonzero only for 
V j^O- That is, indeed, the triplet part of SCDOS only 
contributes to the current if a spin-dependent quasiparti- 
cle distribution is created in the interlayer. The exchange 
field h is chosen to be not very strong h = 10 A. Such 
a choice is in general agreement with the characteristic 
values of the exchange field in weak ferromagnetic alloys. 
However, the results discussed below qualitatively sur- 
vive for the case of more strong exchange fields. Roughly 
speaking, increasing of the exchange field infiuences the 
results in the same manner as increasing of the F layer 
length dp- 
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FIG. 2. Current components (solid line), jgn (dashed line) 
and jln (dotted line) as functions of voltage V , applied be- 
tween the additional eletrodes. The currents are measured in 
arbitrary units. Panel (a) corresponds to dp = 1, div = 2, for 
panel (b) dp = 1, djv = 1 and for panel (c) dp = 2, djv = 1. 
All lengths are measured in units of ^s. The other parameters 
are the following: h = lOA, e'^s = 0.2, T = 0. 



Panels (a), (b) and (c) correspond to different lengths 
of the TV and F regions forming the interlayer. Below all 
the lengths are expressed in units of the superconduct- 
ing coher ence l ength ^g. The magnetic coherence length 
= S,s v^A/ft, is approximately three times shorter than 
^5- For panels (a) and (b) the ferromagnetic layer is not 
thick {dp — !)• They differ by the length of the normal 
layer: panel (a) corresponds to djv — 2 and for panel 
(b) d^r = 1. As it is expected, upon incresing (In the 
magnitude of all the current components decreases not 
very sharply. The corresponding decay length is consid- 
erably larger than ^i?. On the contrary, increase of dp 
suppresses current components js and jsR exponentially 
with the characteristic decay length S^p- It is natural be- 
cause they flow via the singlet and short-range triplet 
components of the anomalous Green's function. These 
components are composed of the Cooper pairs with op- 
posite spin direction and, correspondingly, decay rapidly 
into the depth of the ferromagnetic region. It is seen 
from the figures that for panels (a) and (b) jsR and jlr 
are of the same order, while js is even larger. This is not 
the case for panel (c), where dp — 2. For this parameter 
range js and jsR are already suppressed. However, for a 
certain voltage range (small enough voltages) jLR is not 
suppressed and the dependence of its magnitude on dp is 



the same as on dN- For larger voltages the value of jLR 
is also suppressed. It is interesting to note that this sup- 
pression takes place for all the panels of Fig.[2]irrespective 
of the F layer length. It is obvious that the insensitiv- 
ity of jLR to the length of the ferromagnetic region is a 
result of the fact that it is carried by Cooper pairs com- 
posed of the electrons with parallel spins. However, the 
characteristic behavior of this component in dependence 
on V (sharp maximum at small voltages and subsequent 
suppression) requires an additional explanation. Such an 
explanation is closely connected to the particular shape 
of the anomalous Green's function LRTC and is given 
below upon discussing the LRTC. 




dp/is 
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FIG. 3. Current components (solid line), jg^ (dashed line) 
and jln (dotted line) in dependence on dp/^s (logarithmic 
scale). Panel (a) corresponds to V = 0.05 A, V = 0.1 A for 
panel (b) and V — 0.5 A for panel (c). The other parameters 
are the following: h = lOA, Q'^s = 0.2, djv = 1, T = 0. 

Further, the dependence of all three current compo- 
nents on the length of the ferromagnetic layer is stud- 
ied in more detail. Panels (a), (b) and (c) of Fig. [3] 
demonstrate this dependence for three different voltages 
V. For panel (a) the particular value of this voltage is 
chosen to be = 0.05A. This value approximately cor- 
responds to the maximum of jLR in Fig.[2l For panel (b) 
V = O.lA. Current jLR gradually declines at this voltage 
region. Finally, the plots shown in panel (c) correspond 
to V — 0.5 A, where jLR is already greatly suppressed. 
First of all, it is worth to note that the decay length of 
js and jsR is S.p to a good accuracy for any voltage re- 
gion. Also, it is seen from Fig.|3]that js and jsR oscillate 
upon increasing dp with the period 2tt^p (irrespective of 
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the particular voltage). For js, which is non-zero even for 
spin-independent quasiparticle distribution, these oscilla- 
tions are well-studied. They are a hallmark of mesoscopic 
LOFF-state, as it was mentioned in the introduction. jsB. 
is absent for spin-independent quasiparticle distribution, 
but is carried by the same pairs of electrons with opposite 
spin directions, just as js, and, consequently, also man- 
ifests the LOFF-state oscillations. While the oscillation 
period is the same for js and jsR, there is a phase shift 
between their oscillations, which depends on the partic- 
ular value of the voltage V. 

Unlike js and jsR, jhR does not manifest oscillating 
behavior. Its decay length is not connected to £^f and 
crucially depends on V . This decay length Ilr is maxi- 
mal for the voltage region, where jLR has maximal value 
(Ilr ~ sa for panel (a)) and declines upon in- 
creasing V. 
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FIG. 4. Combinations Fs [panel(a)], Fsr [panel (b)] and Flr 
[panel (c)] as functions of quasiparticle energy e/A. In each 
panel (dp ~ 1, djv ~ 2) for solid curves, (cIf ~ 1, djv ~ 1) for 
dashed curves and {dp = 2, djv — 1) for dotted curves. The 
other parameters are the same as in Fig. [51 

As it was already mentioned in the introduction, the 
dependence of the anomalous Green's function in the in- 
terlayer on the quasiparticle energy can be partially ex- 
tracted from the Josephson current measurements. It can 
be done due to the fact that voltage V enters the cur- 
rent expression (jI9|) only via distribution function (1261) . 
Then, according to Eqs. (HH), ^ and dH]), by taking 
the derivatives of the Josephson currents js, jsR and jLR 



with respect to the voltage applied between the addi- 
tional electrodes, at T — 7> one obtaines that 



djs I dV ^ Im 
djsR/dV ^ Im 
djLR/dV Im 



fliv) 


yA2-(eF4 




fsRiV) 




yA2-(ey4 




fLRiv) 




yA2-(ey + 





(28) 



Here /g^, flj^ and /J are determined by Eqs. ^ 
and ([25]) . respectively. That is, indeed, imaginary parts 
of the anomalous Green's function components coming 
from the opposite interface corresponding to all three 
types of superconducting correlations can be extracted 
from the Josephson current measurements. However, un- 
der the condition jeV^j < A it can be done only for subgap 
energies jej < A. It is worth to note here that the deriva- 
tives of current components js , jsR and j^ji with respect 
to V give us the corresponding anomalous Green's func- 
tion components only in the tunnel limit Gt 1- In 
general case these derivatives are proportional to the ap- 
propriate components of SCDOS, which are expressed via 
the anomalous Green's function in a more complicated 
way. 

Panels (a), (b) and (c) of Fig. 2] represent 
combinations Fs = Im flie)/ \/ A2 — (e -|- iSY 

and Flr, 



Fsr 
Im 



Im 



FsRie)/^' 



work of 



(e -t- i(5)2 , calculated in the frame- 

our microscopic model, in dependence on 
quasiparticle energy e measured in units of A. In each 
panel different curves correspond to different lengths 
dN and dp (See caption to Fig. |4|). It is seen that the 
value of the normal region length does not influence 
qualitatively on all three components of the anomalous 
Green's function. As it is expected, Fs and Fsr are 
strongly suppressed upon increasing of dp- On the 
contrary, Flr is only very weakly sensitive to changing 
of dp- It is dominated by the sharp dip at low energies, 
which followed by wider peaks, where Flr changes sign. 
The width Se of the dip is ~ VADQ'. 

It is the characteristic shape of Flr that is responsible 
for jLR behavior in dependence on V, shown in Fig. [5] 
The point is that at low enough temperatures only the 
part of Flr, belonging to energy interval [— jeT^j, jeFj], 
contributes to jLR- Consequently, upon inreasing of V 
jLR grows sharply up to ^ ~ (l/2)-\/ ADQ' and after that 
starts to decline due to the opposite sign contribution of 
the peaks. It appears that the contributions of the dip 
and the peaks mainly compensate each other, what leads 
to strong supression of jLR for large enough V. The dis- 
cussed above dependence of jLR decay length on V is also 
closely connected to the fact that only the part of Flr, 
belonging to energy interval [— jeFj, \eV\], "works" upon 
creating jLR- Indeed, the characteristic decay length of 
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FLR{e) d'j^{e) |Af(e)| ^ Therefore, Jlr decay length 
~ 1/0' for small voltages and gets shorter for larger volt- 
ages due to increased contribution of higher energies. 
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FIG. 5. Joint currents -I- jg^ (solid curves) and jg + jln 
(dashed curves) as functions of eF/A. For panel (a) dp ~ 
1, djv = 2 and for panel (b) dp — 2, d]v — 1. The other 
parameters are the same as in Fig. O 



IV. S/N/S JUNCTION WITH MAGNETIC 
INTERFACES 



In this section Josephson current is studied for S/N/S 
junction with magnetic S/N interfaces under the condi- 
tion of spin- dependent quasiparticle distribution in the 
interlayer. The model is already described in Sec. |lTl 
The anomalous Green's function in the interlayer is found 
up to the first order in S/N conductance Gt according 
to Eqs. and (IT5|) (assuming h = 0). In general, 
the condensate penetrating into the interlayer region is 
comprised of two types of electron pairs: with opposite 
electron spins and with parallel electron spins. However, 
due to the absence of ferromagnetic elements in the in- 
terlayer region they have the same characteristic decay 
length. The both types of pairs occur in the system if 
there is a reason for spin-flip there. For example, it is 
the case if the magnetization vectors of the both inter- 
faces are not parallel m; H m^. If mi\\mr, then only 
the pairs with opposite electron spins, generated by the 
singlet superconductor, occur in the interlayer region. In 
order to make the formulae less cumbersome we give fi- 
nal expressions only for the case mi\\mr = m. Then, at 
the left {a = +1) and right {a = —1) S/N interfaces the 
singlet part of the anomalous Green's function takes the 
following form 



= fsie' 



/s2e' 



The possibility to extract singlet and triplet compo- 
nents of the proximity-induced anomalous Green's func- 
tion in the interlayer is not the only motivation to study 
the Josephson current under spin-dependent quasiparti- 
cle distribution. Being an easily controllable parameter, 
voltage V gives a possibility to obtain highly nonlinear 
characteristics j{V) with a number of O-tt transitions, 
which can be essential for superconducting electronics. 
As it was already mentioned above, by choosing the ap- 
propriate orientation of Pi and Pr magnetizations, one 
can, in principle, "turn off' either jsR or j^jj current 
contribution. When the full current through the junc- 
tion is given by joint contribution of Jlr and js or jsR 
and js , respectively. The corresponding full currents are 
demonstrated in Fig. [5l Panel (a) reperesents the case 
of short enough ferromagnetic layer dp = ^, while panel 
(b) corresponds to dp = 2. It is seen from panel (b) that 
in this case the main contribution to the current is given 
by Jlr, at least for small enough voltages. It is worth to 
note that it may be not easy to adjust Pi and Pr magne- 
tizations in such a way that only one of the components 
jsR. and Jlr. fiows. In fact, it is not necessary in order 
to obtain highly nonlinear j{V) characteristics. For this 
purpose it is enough to create any spin-dependent quasi- 
particle distribution in the interlayer region. The inter- 
nal structure of the anomalous Green's function Flu can 
also be studied separately for long enough ferromagnetic 
interlayers. 



/si — 



/s2 — 



2iTTGT sinh ef 
ctnZ 

AityGt sinhef 



Aat \cfn 



Z = 4A^ sinh^[AArrfAr] 



Aat — 



1 /G, 



Aw \(7n 



sinh[2AAr(iAr] , 

sinh[Ajvrf7v] , 



G^ 

UN 



(cosh^[AArdAr] + l) + 



A \9a\^ sinh^ [AAfrfAf] 
' ' A? 



(29) 



where Ajv is determined below Eq. (|A4p . 

The triplet component of the anomalous Green's func- 
tion has only z-component and takes the form 



//' = (0,0,/,) 



fz = hie 



-iax/2 



iax/2 



47rG'rf,GT sinh 65 
fzl = ^^^r^ ^ X 



1 /G, 



Aat \<yN 



sinh^[AAr(iAr] + 2 



S^G^Grsinhef , , 

fz2 = ^ cosh[A7va7vJ 



(30) 
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where Z is determined in Eq. ((29|) . Physically, fs,z\ are 
generated by the proximity effect at the same S/N in- 
terface and js.z2 are extended from the opposite S/N 
interface. 

Just as in the previous section, in order to generate 
a spin-dependent quasiparticle distribution in the inter- 
layer, additional electrodes are attached to it. The prin- 
cipal scheme is the same as before except for the fact that 
there is only one normal region in the considered system. 
Therefore, we assume that the interlayer is attached to 
two additional normal electrodes and Nt and electrode 
N-h has insertion P made of a strongly ferromagnetic ma- 
terial. The unit vector aligned with the magnetization 
of P is denoted by M. Again, if voltage 2V is applied 
between the electrodes iVj, and iVt, then the electric po- 
tentials for spin-up and spin-down electrons in the N-^ 
region, inclosed between P and the normal interlayer, 
counted from the level of the superconducting leads are 
n = - 14)/2 = V and n = {^t ~ Vb)/2 = -V . The 
distribution functions for spin-up and spin-down elec- 
trons in this region are close to the equilibrium form 
(with different electrochemical potentials). In a matrix 
form the distribution function is expressed by Eq. (j26p 
with the substitution M for Mi. 

Now we can obtain the distribution function in the 
interlayer, which enters the current [Eq. (|19p ]. Again, 
for simplicity we assume that <C 1. Consequently, 
the dissipative current flowing through Nt/N/Nt junc- 
tion is negligible and, therefore, the y-dependence of the 
distribution function in the interlayer region can be dis- 
regarded. For simplicity we assume below that 7Vf||m. 
Under this condition the distribution function i^*^*^-* in the 
interlayer calculated according to Eq. (ppj) at h = sup- 
plemented by boundary conditions at S/N interfaces ((2T|) 
is spatially constant and equal to its value coming from 
iVb region. If M i/( m, then the spatially constant dis- 
tribution function does not satisfy boundary conditions 
(1211) any more. In this case the problem become two- 
dimensional and much more complicated. 

Now we are able to calculate the Josephson current 
through the junction according to Eq. ([19]). After substi- 
tution of the expression for the singlet part of the anoma- 
lous Green's function [Eq. d^H]) ] and the scalar part of the 
distribution function [Eqs. ([25]) and (fT5)) ] into first two 
terms of Eq. (fT9)) . the contribution of the SCDOS singlet 
part takes the form 



stitute Eqs. dSO]), ([26]) and (HSl) into Eq. 



Jt 



iC^G^ sinx 



ea 



N 



A'^deiptis) X 
cosh[AjvdAr] 



(32) 



As opposed to the problem of S/NFN/S junction con- 
sidered in the previous section, it is seen from Eq. (j32p 
that jt values at the left and right S/N interfaces are 
equal to each other. As for the case of S/NFN/S 
junction, the part of current ([T9l) generated by the 



term cx cosh <d§ 



(0)/ 



/2 vanishes due to 



the fact that the scalar part (p^ of the distribution 
function in the interlayer [Eq. (|26p ] is an odd func- 
tion of quasiparticle energy. Further, under the condi- 
tions \eV\ < A and T <C A the last term, generated 
by cx aGA/flCOshefm [(^("'(e) -t- (p'-°^-e)] /2 also van- 
ishes because this expression is an odd function of quasi- 
particle energy at |e| < A and is absent elsewhere. Tak- 
ing into account that m/||mr||iVf one can obtain from 
Eq. dH]) that d^ip'^^^ = at the S/N interfaces. There- 
fore, (^(^) is approximately constant in the interlayer. 
Moreover, this constant is to be equal to zero in order to 
satisfy the condition — j^. Therefore, the full Joseph- 
son current flowing through the junction is given by the 
sum of singlet [Eq. §1^] and triplet [Eq. SCDOS 
contributions. 
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(b) 



2iG^ sin Y / . 9 , , , , 
3s = / A^depoie) x 



ecrjv 



Ajv sinh[AjvdAr](l 



[(£ + Z(5)2 - A2]Z(£) 



(31) 



where Z{e) is determined in Eq. ([29)) . 

The current flowing through the SCDOS triplet part 
and expressed by the third term in Eq. ([T9|) takes the 
form [in order to obtain this expression one should sub- 



FIG. 6. (a) Functions Fs (dotted line) and Ft (solid line) as 
functions of e/A for S/N/S junction with magnetic interfaces, 
(b) Full critical current (solid line) and its contributions j'J 
(dotted line) and jf (dashed line) as functions of eV/A. For 
the both panels djv = 0.5^s, Gtf.^s/o'N = 0.35 and T = 0.1 A. 

As for the previous case of S/NFN/S junction, is an 
even function of voltage V applied to the additional elec- 
trodes and jt is an odd function of this voltage. There- 
fore, contributions j^ and jt can be extracted from an 
experimentally measurable Josephson current and, so, it 
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makes sense to discuss them separately. Panel (b) of 
Fig. ini demonstrates the full critical Josephson current 
and its contributions jg and as functions of V for a typ- 
ical set of parameters (See caption to Fig. ^ for specific 

values). Functions Fs{e) 



Im 



and Ft{e) = Im 



- (e + iS)^ 
are repre- 



sented in panel (a) of Fig. |6] for the same set of parame- 
ters. As it is seen from the definitions given in Eqs. 
and (|30l) . these functions are proportional to the singlet 
and triplet components of the anomalous Green's func- 
tion, coming from the opposite S/N interface, and can 
be experimentally found by differentiating the currents 
and with respect to voltage as it was explained 
in the previous section. 
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FIG. 7. Fs and Ft as functions of e/A for S/N/S junc- 
tion with magnetic interfaces. The upper row represents Fs, 
while the lower row demonstrates Ft. For panels (a) and (c) 
djv = different curves correspond to different values 

of G^^s/crjv — 0.1 (black solid curve), 0.3 (dotted curve), 0.7 
(dashed curve) and 1.1 (gray solid curve). In panels (b) and 
(d) G^(,s/o'N = 0.5 and different curves correspond to dif- 
ferent d]v/Cs ~ 2 (black solid curve), 1 (dotted curve), 0.6 
(dashed curve) and 0.4 (gray solid curve). T = 0.1 A. 

The characteristic shape of Fg and Ft dictates how js 
and jt behave in dependence on V. Upon discussing the 
characteristic features of js and jt we only consider V > 
and, correspondingly, e > for F^ and Ft. The main 
characteristic features of Fg and Ft , which are responsible 
for the current behavior, are proximity induced dips at 
£0 ^ GffyS^gA/aNdN (for the parameter region < A). 
These dips are followed by abrupt changing of sign of the 
corresponding quantity. Fig.[7]shows F^ and Ft evolution 
with 6*0 (left column) and with (right column). It is 
seen that upon increasing the proximity induced dip 
shifts to higher energies. If the junction becomes shorter 
the dip also shifts to the right and its integral height 
increases due to the fact that the proximity effect is more 
pronounced for short junctions. 

According to Eqs. and (pS)) at low enough temper- 
atures only the part of Fg, belonging to energy intervals 
[— oo,— |ey|] and [|ey|, -f cx)], contributes to j^. Conse- 
quently, upon inreasing of V the absolute value of 
grows up to y ~ £0 and after that starts to decline due 



to the sign changing of Fs(e) at e = e^. Analogously, only 
the part of Ft, belonging to energy interval [— |ey|, \eV\], 
contributes to jt- Therefore, the absolute value of jt 
also grows up to V ^ and declines after that. The de- 
scribed behavior is characteristic for the absolute value of 
j^ and j^ as for eV > 0, so as for eV < 0. However, due 
to the fact that j^ is symmetric and jf is antisymmetric 
function of V, the total Josephson current is highly non- 
symmetric with respect to as it is seen in Figs. |6] and [H 
While for eV < the contributions of and jt partially 
compensate each other leading to suppression of the full 
current and — 7r-transition at some finite V, they are 
added for eV > resulting in the considerable current 
enhancement. The value of eV, where the peak in the 
critical current is located, can be used for experimental 
estimate of spin-mixing parameter G^, characterizing the 
magnetic interface, because eT^ ~ e^. For short enough 
junctions with < ^5 the current value can even ex- 
ceed the critical current value for S/N/S junction with 
nonmagnetic S/N interfaces (G^, — 0) and the same S/N 
interface conductance Gt for some voltage range. It is 
worth to note here that such an enhancement is only pos- 
sible for finite V, when the triplet part of SCDOS Ft con- 
tributes to the current. At V = the critical Josephson 
current through S/N/S junction with magnetic interfaces 
Gip ^ is always lower than the corresponding current 
fot S/N/S junction with nonmagnetic intefaces but the 
same interface conductance Gt (this statement is valid 
for the entire range of parameters we consider) . 
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FIG. 8. Full critical Josephson current for S/N/S junction 
with magnetic interfaces in dependence on eV/A. Panel (a) 
demonstrates the case of low-temperature junction, where 
the proximity effect is well-pronounced: T — O.OIA, djv = 
O.S^s. G^.^s/o'N = 0.15(dashed curve), 0.3 (solid curve) 
and 0.45 (dotted curve). Panel (b) corresponds to longer 
junction at the same temperature: T — O.OIA, cIm = S^s- 
G^i^s/o'N ~ 1 (dashed curve), 2 (dotted curve), 3 (solid curve) 
and 4 (dashed-dotted curve). Panels (c) and (d) represent the 
same results as panels (a) and (b), respectively, but at higher 
temperature T = 0.1 A. For all the panels the gray solid line 
represents jnm{V) for the corresponding set of the parame- 
ters. 

Let us denote the value of the critical current for S/N/S 
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junction with nonmagnetic intefaces and interface con- 
ductance Gt by jnm{y)- In the framework of the mi- 
croscopic model of S/N interface considered in Appendix 
IbI comparison between j{V) and inm.{V) physically cor- 
responds to comparison between the Josephson currents 
in the system with a thin magnetic layer between N and 
I and without it. The value jnm{V = 0) is shown in 
Fig. IHJb) by the horizontal line. It is seen that a small 
excess of the total current j'^ over jnm{V = 0) takes 
place for some voltage range. However, the excess can 
be much greater, the current at finite V can exceed the 
equilibrium current jnmiV = 0) for nonmagnetic S/N/S 
junction more than twice. Such a case is illustrated in 
Fig.mja). Maximal excess can be expected for short junc- 
tions with £0 « A, where the proximity effect in Ft is 
most pronounced and the proximity induced dip at 
merges with the coherence peak at A thus greatly en- 
hancing Ft value in the subgap region. In addition the 
temperature should be low enough in order to avoid tem- 
perature smearing of the effect. 

For longer junctions with cIn ^ S,s the current does 
not exceed the equilibrium value jnm{V — 0) because 
of weaker proximity effect in the interlayer region, as it 
is ihustrated in Fig. [DJb). For all the panels of Fig. [8] 
the gray solid line represents jnmiy) for the correspond- 
ing set of the parameters. It is worth to note here that 
dependencies inmiV) on V are qualitatively very simi- 
lar to the current discussed in Ref. 30 for nonmagnetic 
S/N/S junction under nonequilibrium quasiparticle dis- 
tribution in the normal interlayer. Indeed, at — 
triplet part of SCDOS Ft is absent and, consequently, 
the vector part of the distribution function (l26l) does not 
contribute to the current. The singlet part of this distri- 
bution function is formally equivalent to the nonequilib- 
rium distributioEk^ for a narrow normal interlayer (a wire 
or a constriction). Full quantitative agreement between 
our results for jnm{V) and the results of Ref. [13 can- 
not be reached because they are obtained for somewhat 
different parameter ranges. 

Panels (c) and (d) of Fig. |S] show the results for the 
current at higher temperature T — 0.1 A. Panel (c) 
corresponds to shorter junction with djsi = 0.3^5, while 
panel (d) demonstrates the case of longer junction with 
disi = 3^s. It is seen that for short junction, where 
the effect of current enhancement is well pronounced at 
low temperatures, raising of the temperature suppresses 
the effect. The reason is that distribution function (p6)) 
smears upon raising of the temperature. Consequently, 
not only the part of Ft corresponding to |e| < e^, but 
also some region of higher energies, where Ft has oppo- 
site sign, is involved into the current jt now. This leads 
to partial compensation of Ft parts with different signs. 

Although the maximal value of the critical Josephson 
current, which can be reached at a finite V , is suppressed 
by temperature, the current dependence on T at a par- 
ticular voltage V can be quite interesting. Fig. 1^] demon- 
strates how the current depends on temperature at sev- 
eral specified values of voltage V for the case of short 




FIG. 9. Full critical Josephson current for S/N/S junction 
with magnetic interfaces in dependence on temperature for 
several different voltages V. = O.S^s, black solid curve: 
G^£,s/<yN ~ 0.3, eF/A = 0.85; dashed curve: G^£,s/<jn = 
0.3, eV/A = 0.7; dotted curve: G^is/(yN = 0.3, eV^/A = 
0.5; dashed-dotted curve: G^is/oN = 0.15, eV/A = 0.83. 
Gray solid curve represents the temperature dependence of 
jnmiy ~ 0) for the corresponding set of parameters. 



junction with djsi = 0.3^5. The gray solid line represents 
the dependence of jnmiV = 0) on temperature and is 
given for comparison of our results with the equilibrium 
nonmagnetic case. It is well known, that the Josephson 
current for equilibrium nonmagnetic S/N/S junction de- 
clines upon raising of temperature rather sharply, as it 
is demonstrated by the gray solid curve. On the con- 
trary, the current at finite V for a S/N/S with magnetic 
interfaces can even grow up to some temperature and 
only after that start to decline. The qualitative expla- 
nation of this fact is the following. The main contribu- 
tion to the Josephson current in nonmagnetic equilibrium 
S/N/S junction is given by high peak of SCDOS located 
at low energies. Consequently, the temperature smearing 
of the equilibrium distribution function tanh£/2T cru- 
cially reduces the current. At the same time, the main 
contribution to jt is given by the energies up to e ~ \eV\ 
and under the condition that \eV\ < the temperature 
smearing of the distribution function involves higher en- 
ergies, where the absolute value of Ft even larger, in the 
current transfer. In addition, at some voltage ranges the 
junction can manifest — tt transition in dependence on 
temperature. 



V. SUMMARY 

In conclusion, we have theoretically investigated the 
Josephson current in weak links, containing ferromag- 
netic elements, under the condition that the quasiparti- 
cle distribution in the weak link region is spin-dependent. 
Two types of weak link are considered. The first system 
is a S/N/F/N/S junction with complex interlayer com- 
posed of two normal metal regions and a middle layer 
made of a spiral ferromagnet, sandwiched between them. 
The second considered system is a S/N/S junction with 
magnetic S/N interfaces. In both cases spin-dependent 
quasiparticle distribution in the interlayer region is pro- 
posed to be created by attachment of additional elec- 
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trodes with ferromagnetic elements to the interlayer re- 
gion and applying a voltage V between them. Interplay 
of the triplet superconducting correlations, induced in 
the interlayer by the proximity with the superconduct- 
ing leads, and spin-dependent quasiparticle distribution 
results in the appearence of the additional contribution 
to the Josephson current jt , carried by the triplet part of 
SCDOS. 

It is shown that jt is an odd function of V, while the 
standard contribution j^, carried by the singlet part of 
SCDOS, is an even function. So, jt can be extracted 
from the full Josephson current measured as a func- 
tion of V. Further, it is demonstrated that deriva- 
tive djt/dV can provide direct information about the 
anomalous Green's function describing the superconduct- 
ing triplet correlations induced in the interlayer. We 
show that in the S/N/F/N/S junction the contributions 
given by the short-range (SRTC) and long-range (LRTC) 
components of triplet superconducting correlations in the 
interlayer can be measured separately. 

For S/N/S junction with magnetic interfaces it is also 
obtained that the critical Josephson current at some 
finite V can considerably exceed the current flowing 
through the equilibrium nonmagnetic S/N/S junction 
with the same S/N interface transparency. This enhance- 
ment is due to the fact that the triplet component of SC- 
DOS "works" under spin-dependent quasiparticle distri- 
bution giving the additional contribution to the current, 
while it does not take part in the current transfer for spin- 
independent quasiparticle distribution. In addition, we 
have studied temperature dependence of the critical cur- 
rent in S/N/S junction with magnetic S/N interfaces. As 
opposed to the case of equilibrium nonmagnetic S/N/S 
junction, where the current is monotonouosly suppressed 
by temperature, in the considered case at a finite volt- 
age V it can at first rise with temperature and only then 
start to decline. 

The dependence of the full critical current on V is typ- 
ically highly nonlinear and strongly nonsymmetric with 
respect to = due to the interplay of js and jt- This 
also leads to appearence of a number of O-tt transitions 
in the system upon varying controlling voltage V . 



Appendix A: Microscopic calculation of anomalous 
Green's function and Josephson current in 
S/NFN/S junction 

In this Appendix we calculate the anomalous Green's 
functions fs, fsR and /lr and the corresponding cur- 
rent contributions js, jsR and j^B, in the framework of 
the most simple microscopic model for N/F/N interlayer. 
We assume the NF interfaces to be absolutely transpar- 
ent. This approximation simplifies the calculations sig- 
nificantly, but does not influence qualitatively our main 
conclusions. For this case the boundary conditions at 
X = take the form 



(^NdxQN — crpdxgF 



(Al) 



As far as we only need anomalous Green's functions to 
the first order in S/N interface transparency, the above 
boundary conditions should be linearized. Then for re- 
tarded and advanced Green's functions they read as fol- 
lows 



fR,A _ fR,A 
In ~ JF 



IN 



OfOxJp 



(A2) 



The boundary conditions for the distribution function 
at N/F interface to the considered accuracy take the form 



(fF = (fiN , 

ardxipF = cTNdxipN 



(A3) 



The singlet part of the anomalous Green's function, 
calculated according to Eqs. (|13p . ([T5|) and (jA2p . at the 
left {a = +1) and the right (a = —1) S/N interfaces 
takes the following form 



tanh (pN sinh O^e 



cos(x/2) cosh(/)+ 
A+ sinh (/)_|_ -|- p cosh (/)_|_ 

cos(x/2) cosh0_ 
A- sinh ( 



R^-^ax/2 ^ ^TtCt sinh Of 

2a F cosh^ (t>N 
ia sin(x/2) sinh 0+ 
A+ cosh 0+ + p sinh 0+ 
ia sin(x/2) sinh 0_ 
p cosh A_ cosh 0_ + p sinh 



(A4) 
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where A± = ^/h/D^l^i), Xn = \/-2i{e + iS)/D, (j)± = 
X±dF/2, 4>N = ^Nd^/^ and p = {a^ /o'F)XNtanh4>N. 

The results for triplet components fsR and Jlr are 
the following 
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A-i- sinh + p cosh ( 
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A_ sinh (/)_ + p cosh 



(A5) 



where At = y ^ ^ 2i(e + iS)lD and (jit = \tdF/2. 

The fact that /s^ rapidly decays in the ferromagnetic 
region and, consequently, represents the SRTC can be 
easily seen from Eq. (|A5|) in the limit of thick enough F 
layer: dp/^F ^ 1- To the leading order in the parameter 
g-^F/CF foj- quantities /^^ and /g^, defined by Eq. 
one obtains from Eq. (jASP 
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fsR. — ^' 



IttGt sinh 
aF cosh^ (j)N 



R r 



A 



+ e 



A+ + P 



A- 



(A++p)2 (A_+p)2 



(A6) 



At the same regime the corresponding components of 
/lr. take the following form 



I 

LR 



R r 



ittGt sinh 0^ 
2a F cosh^ (f>N 
e'(At cosh[2(/.t]- 



1 



1 



A++P A_ 
p sinh [2 



(p cosh (/>t + At sinh 4>t)ip sinh 0t + At cosh (j)t ) 



fr 

JLR ~ 



inGr sinh Qg 
2a F cosh^ 



1 



1 



On 

e'At 



A++P A_ 



(pcosh^t + At sinh 0t)(/3 sinh 0t + At cosh(/)t) 



.(A7) 



As it is seen, does not contain the small factor 

g-dp/^F ^YiQ leading approximation and, therefore, Jlr 
describes the LRTC. The characteristic decay length of 
fLR in the F layer is |At|"^ 

To the considered accuracy the singlet component of 
the anomalous Green's function also decays at the dis- 
tance ~ £_F in the F layer, just as the SRTC fsR does. In 
the regime dp/^F ^ 1 it can be obtained from Eq. (|A4|) 
that 



iirGi 



iirGx sinh Qg ^ 



tanh 0JV sinh Qg 



2a F cosh (/)jv 
inGr sinhO^ 



1 



1 



ap cosh 



6 AT 



A+ + p A_ + p_ 



A+e 



(A++p)2 (A-+p)2 



.(A8) 



Substituting Eqs. (jA4p and (jASp together with the ex- 
pressions for the vector part of the distribution function 



[Eqs. (Uni) and ([T8|)] into Eq. ^ one can find 

oo 

. _G|,sinx [ iA'^de(fo{e) 



Seap J [{e + iS)^ - A2] cosh^ 0jv 



1 



tanh ( 



A_|_ tanh 0+ + pn ^+ + Pn tanh 0+ 
1 tanh (t>_ 



JSR 



A_ tanh 0_+ Pat A_ + pat tanh 

oo 

G^sinx f iA'^dsiptis) 



(A9) 



Seap J [A2 - (e + iS)^] cosh^ 



9 AT 



1 



tanh ( 



A_|_ tanh 4)+ + pn ^+ + Pn tanh 
1 tanh 



A- tanh (f>^ + pN A_ + pA? tanh 

oo 

Gf^sinx f iA^de(pt{e) 



, (AlO) 



8ecr_F 7 [A2 - (e + cosh^ 0Ar 



6' 



p -I- At tanh (/)t 
tanh fli- 



tanh ( 



A_ + ptanh( 
1 



A+ + p tanh 0+ 
0' tanh 0t 
p tanh 0f -t- At 
1 



A+ tanh 0+ + p A_ tanh 0_ -f p 



(All) 



Appendix B: Microscopic model of magnetic S/N 
interface 

Let us introduce the electronic scattering matrix S'^ 
associated to electrons with spin a of the n*^ transmission 
channel. We assume that the interface do not rotate an 
electron spin, that is the scattering matrix is diagonal in 
spin space. 



ncr ncr 
+/ r'' 



(Bl) 



where Tna denotes the reflection amplitude at the left 
(right) side of the interface and tn^'' the transmission 
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amplitude from the left (right) side to the right (left) 
side of the interface. Taking into account the constraints 
on S"^ resulting from the unitarity condition S'^S'^^ — 1 
and time reversal symmetry one can show that without 
any loss of generality S'^ is entirely determined by the 
following parameters: the transmission probability T!„, 
the degree of spin polarization P„ and the spin-mixing 
angle d(pn^\ These parameters are defined as T„o- = 
\tna\^ = r„(l + (7P„) and arg[r,L^;^] = <ff'>+a{d^f^/2). 

These parameters can be straightforwardly calculated 
in the framework of a microscopic model describing the 
interface. Here we model S/N interface by an unsulating 
barrier I (with a transparency T 1) and a thin layer 
of a ferromagnetic metal, which is located between I and 
the normal inter layer. This layer is supposed to provide 
a required value for the spin-mixing angle. Given the 
exchange field in the ferromagnetic layer is small with 
respect to the Fermi energy h <^ ep, in the framework 
of this microscopic model Tn = T , Pn ^ Q and dtpn ~ 
2wFh/vF, where wp '^s the length of the ferromagnetic 
layer and is the corresponding Fermi velocity. 

The main parameters entering magnetic boundary con- 
ditions ([9]) are connected to the microscopic parameters 
r„, Pn and dip'n^'^ by the following way^ 

Gt5 = 2G,^T„ , (B2) 



GmrS — Gq'^^TnPn , (B3) 



G^S = 2G, ^(T„ - l)dipn , (B4) 

n 

where S is the junction area and Gq = /h is the quan- 
tum conductance. It is worth to note here that bound- 
ary conditions (jH]) are the expansion in small T„ , P„ and 
dipn of the more general boundary conditions'*'*' and, con- 
sequently, Eq. ([9]) is only valid if all these parameters 
are considerably less than unity. However, because of 
summation over large number of transmisson channels, 
it does not mean that the parameters Gt, Gmr and G^ 
must be small. Let us estimate the value of dimensionless 
G^ = G^^s/<^N, which can be obtained in the framework 
of our microscopic model. For T <^ 1 



Sa^ T ' ^ ' 

where N is the number of transmission channels and / is 
the mean free path, dip means the average value of the 
spin- mixing angle c?(p„ . For rough estimates it is possible 
to take dip ~ 2w ph / v f ■ 

Our main results for S/N/S junction with magnetic 
interfaces are calculated for G^ ^ 1. From Eq. (IBSP it is 
seen that it is quite reasonable to expect such values of 
for a magnetic interface composed of an unsulating 
barrier and a weak ferromagnetic alloy with wp ^ Cf- 
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